Spreading of thin films assisted by thermal fluctuations 
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We study the spreading of viscous drops on a solid substrate, taking into account the effects 
of thermal fluctuations in the fluid momentum. A nonlinear stochastic lubrication equation is 
derived, and studied using numerical simulations and scaling analysis. We show that asymptotically 
spreading drops admit self-similar shapes, whose average radii can increase at rates much faster 
than these predicted by Tanner's law. We discuss the physical realizability of our results for thin 
molecular and complex fluid films, and predict that such phenomenon can in principal be observed 
in various flow geometries. 
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Water drops spreading on a table and oil drops lubri- 
cating a pan are two common examples of a phenomenon 
encountered frequently in the kitchen as well as in nat- 
ural and industrial environments: Spreading of liquids 
on solid surfaces. Despite its prevalence and the basic 
hydrodynamic principles involved, it was not until the 
late 1970's that the asymptotic rate of spreading pro- 
cesses was found by Tanner 1] for surface-tension domi- 
nated flows. The spatial scale t of a viscous drop spread- 
ing on a smooth plane increases asymptotically in time 
as I ~ t z , where z — 1/10,1/7 for radially symmet- 
ric two-dimensional and one-dimensional flow geometries, 
respectively. This asymptotic response has been found 
in many molecular and polymeric drops, whose decreas- 
ing thickness has been detected down to 100 nm 0- 
As the field of nano-fluidics is evolving towards forma- 
tion of thinner and thinner liquid films, theoretical tools 
are needed to describe flow patterns in such geometries. 
However, the applicability of classical hydrodynamic the- 
ory for these systems is questionable. While the necessity 
to incorporate van der Waals (vdW) fluid-solid attraction 
was recognized long ago 0, other fundamental aspects 
have never been fully addressed. In particular: Does a 
three-dimensional hydrodynamic description hold for a 
film whose thickness is just few molecular layers? What 
are the effects of thermal fluctuations on the determinis- 
tic hydrodynamics at such small scales? 

To resolve these questions, extensive molecular dynam- 
ics (MD) or lattice-gas (LG) simulations of flow in liquid 
films are required to allow comparison and quantify de- 
viations from a regular hydrodynamic theory While 
a full resolution of these questions is still not available 
within current computational possibilities, a few recent 
studies are indicative. For example, Abraham et al. |5j 
used LG algorithm to study flow in a precursor film as- 
sociated with a spreading drop, and demonstrated signif- 
icant deviations from the predictions of a hydrodynamic 
model On the other hand, MD simulations of nano- 
jets, whose initial radius was around 10 molecular di- 
ameters, were shown to be qualitatively consistent with 
simulations of a stochastic Navier-Stokes (NS) equation, 



where the viscous stress tensor was supplemented by a 
stochastic tensor whose temperature-dependent magni- 
tude is determined from the fluctuation-dissipation the- 
orem |f|. The emerging picture is that, at least in some 
cases, a hydrodynamic description can still be used as a 
quantitative tool in studying flow of nano-fluids, however 
modifications of the classical equations are required. 

In this Letter we take one further step forward by ex- 
ploring the influence of thermal fluctuations on the shape 
and rate of spreading of nano-dimension drops, while 
assuming a generalized NS equation holds, similarly to 
0. We should note, however, that in addition to ther- 
mal fluctuations other modifications of the hydrodynamic 
equations that stem from density variations near the in- 
terfaces 0, might be necessary in this regime. 

Let us start by considering the dynamics of the height 
h(x,y,t) of an infinite incompressible planar viscous fluid 
film on top of a smooth solid surface, located at h = 0, 
as depicted in Fig 1. We consider highly viscous fluids, 
such that inertia can be neglected. The mass conserving 
dynamics of long wave length fluctuations of the surface, 
\Vh\ <C 1, is described by the lubrication equation 
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where rj is the viscosity and p is the pressure. The deriva- 
tion of Eq. from the NS equation is a standard ex- 
ercise in fluid mechanics ||. Spatial variations of the 
pressure associated with fluctuations of the liquid-vapor 
interface result from several sources: gravity, surface ten- 
sion, and VdW attraction with the solid surface: 



p = pgh~ jV 2 h + A/h 3 , 



(2) 



where p is the fluid density, g is the gravitational ac- 
celeration, 7 is the liquid-vapor surface tension, and A 
is the Hamaker constant. Tanner's law corresponds to 
similarity solutions of Eq. Q 



h(S,t) = \x\-PfQx\/t*) 
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in the surface-tension-dominated regime |7V 2 ft.| ^> 
\pgh\, \A/h 3 \, where the exponent (3 = 1, 2 is determined 



2 



by requiring volume conservation V — J d d xh(x,t), 
yielding /3 = 1 for one-dimensional (d — 1) and (3 = 2 for 
two-dimensional (d — 2) flow geometries, respectively. 

The central equation of this paper is a stochastic gen- 
eralization of Eq. Q : 



dh 1 .,„_ . /2fc B T, 



V-[/i 3 ^(f,t)] , (4) 



which captures effects of thermal fluctuations on the sur- 
face dynamics. Here, £(a;, t) is a spatio-temporal Gaus- 
sian white noise. Principally, Eq. 10} can be derived 
from the full three-dimensional NS equation, similarly to 
Eq. by adding a stochastic stress, representing ther- 
mal fluctuations of the fluid momentum, to the viscous 
stress tensor [fj- We can avoid however such a tedious 
derivation by considering a reduced, linear version of Eq. 
(|T|). and use the fluctuation-dissipation theorem to find 
the correct magnitude of a Langevin term that gives rise 
to equipartition of the thermal energy density carried by 
its eigenmodes h^(x,t) = H + Sh^(t) cos(<f • x). Here 
\8hqj H\ <C 1, H is the average thickness of the film, and 
q is a planar wave vector. Namely, the linear eigenmodes 
of the stochastic surface dynamics are required to satisfy: 
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where Vq = (P9 + 7l?| 2 — 3A/ H 4 )6h^cos(q ■ x) is the pres- 
sure (mechanical energy density) of a surface eigenmode, 
T f j — 3i]/\q\ 2 H 3 is its friction coefficient, and £$(t) is 
the spatial Fourier transform of £ (a?, t) . Notice that the 
long wavelength approximation underlying Eq. im- 
plies \q\H <C 1. Dividing both sides of Eq. (JHJ by T| g | 
and taking the inverse Fourier transform we obtain: 
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The linear Eq. JSJ describes near equilibrium thermal 
fluctuations of a surface, \6h\ <C H. The spreading dy- 
namics of a drop that does not satisfy this condition must 
be described by a nonlinear equation. To this end, no- 
tice that Eq. can be recovered from the deterministic 
part of Eq. © by making the transformation H,Sh — > h 
and requiring the resulting equation to conserve fluid 
mass. By following exactly the same steps the nonlin- 
ear Langevin Eq. 0} can be derived from Eq. ©. 

A word of caution is in order here. As for any ex- 
trapolation of the fluctuation-dissipation theorem to non- 
linear, far from equilibrium dynamics, a local equilib- 
rium assumption must be made It}]. Namely, the de- 
scription of the surface dynamics with Eq. Q assumes 
that the magnitude of thermal fluctuations of the liquid- 
vapor interface at x is determined solely by the local 
value of the surface height h(x,t). This assumption is 
justified only for thermal fluctuations whose wavelength 




FIG. 1: Schematic representation of a one-dimensional 
spreading drop, confined in a channel of width W . The pre- 
cursor layer of height b is also shown. 



A <C A*(x, t) = h/\Vh\. The relaxation dynamics of fluc- 
tuations whose wavelength A > A* (x, t) is strongly cou- 
pled to interface fluctuations, and thus their magnitude 
cannot be assumed to be given by the near-equilibrium 
result (@J. Including these modes in the stochastic anal- 
ysis involves advanced methods [T(il | , and will not be pur- 
sued here. Since we expect \*(x,t) — > oo as t — > oo, 
our local equilibrium assumption becomes asymptotically 
correct. In addition, the number of linear eigenmodes 
with A > A* scales as (A*/L) where L is the lateral 
size of the drop. Therefore, we expect Eq. 0} to provide 
a better description of the pre-asymptotic dynamics in 
two-dimensional than in one-dimensional geometries. 

In studying Eq. 0} , our basic motivation was to under- 
stand the possible effects thermal fluctuations may have 
on the asymptotic rates of spreading, e.g. by modifying 
Tanner's law. With this view, we focused our analysis on 
two characteristic flow geometries: (i) one-dimensional 
drops confined in a channel of width W (Fig. 1), and (ii) 
two-dimensional radially symmetric drops. In the first 
case, Eq. Q assumes the form: 



dh 7 d /^3^ 3 h 
dt 37/ dx dx 3 
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In the second case, Eq. (@J becomes: 



h s ' 2 £(x,tj\ . (7) 
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We simulated the volume-conserving spreading of a drop 
dominated by the dynamics J7} by using finite-difference- 
based computational techniques as in Uil, which guaran- 
tee non-negativity of the field h(x,t) for the deterministic 
part of Eq. J7}. In our simulations, the noise term is in- 
cluded in the RHS of a spatial-temporal discrete version 
of Eq. , which is advanced in time through an implicit 
method. Non-negativity of h(x, t) is physically imple- 
mented by a short-range repulsive potential between the 
liquid and the solid substrate, which in principle should 
be included in the pressure in Eq. (0} . We avoid the ex- 
plicit use of such a potential, by allowing numerical noise 
realizations only if they preserve non-negativity oih{x 1 1). 
Such a procedure induces correlations in the otherwise 
white noise field, which are unavoidable if the repulsive 
potential is not introduced explicitly. In our simulations 



3 




Log(t) 



FIG. 2: Average shape of the droplet (h(x, t)) over 50 real- 
izations as a function of time for zero temperature (left) and 
a = 1CP 2 (right). The right picture clearly shows the assisted 
spreading of the droplet in presence of thermal fluctuations. 
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we use a non-dimensional version of Eq. © : 



9/i 



9 ,r,a 3 /i, 



9 



^/ 2 e(£,t) 



(9) 



where h = h/hg , x — x/ho , t = t/toi where /io is 
the maximal height of the initial drop, t Q = 3r)h /"f, 
and a = ksT/^Who. Our initial condition is a one- 
dimensional droplet with circular cross section, and a 
substrate wet by a precursor film of height b to ensure 
complete wetting (see Fig. 1). Physically, b represents a 
scale below which the surface dynamics is not governed 
by surface tension. For molecular liquids, this scale is 
related to the dominance of VdW forces near solid walls, 
whereas for complex fluids it is determined by the size of 
the fluid constituents a (e.g. for colloidal suspensions a is 
the colloid diameter). The results of our simulations are 
independent of 6, as long as b is small enough compared 
to other length scales involved in the dynamics ©. 

Typical drop shapes obtained by averaging h(x,t) over 
many realizations of Eq. © with several values of a are 
presented in Fig. 2, and compared to the dynamics of 
a spreading drop governed by the deterministic Eq. IjTjI. 
Fig. 2 clearly indicates the effect of thermal fluctuations 
in enhancing the rate of spreading drops. 

To gain some quantitative understanding of this effect, 
we measured the average rate by which characteristic lat- 
eral scales I of the drop evolve for various magnitudes of 
the stochastic force, and compared them to Tanner's law 
in this geometry: Idet ~ t 1 / 7 . To estimate £ we have used 
the averaged second moment of the height profile, 



£ = 



— I dx(x — x) 2 h(x, t) 



-|l/2\ 



(10) 



where x = [J dx xh{x,i)]/V is the instantaneous droplet 
center position, V = J dx h(x, t) is the constant volume 
of the droplet, and (• • • ) represents an average over real- 
izations of the noise £(x,i). The results of this analysis 
are shown in Fig. 3a, from which we extract a modi- 
fied asymptotic spreading rate in volume-preserving one- 
dimensional flow geometry: £ s toch ~ t 1 / 4 . Obviously, the 
larger a is the earlier is the deviation from Tanner's to 



FIG. 3: Results of volume conserving dynamics in one- 
dimensional geometry, a) Lateral scale of the droplet (I1UB 
as a function of time. Solid lines represent averages over 50 
realizations of @ with a = 10 -2 , 5 x 10~ 3 , 10 -3 (from top to 
bottom), while the dashed line is the noiseless (a = 0) dy- 
namics. Dotted lines correspond to the power laws P^ 4 and 
t 117 . b) Averaged profile of the droplet for times t = 10 3n/4 , 
n = 0, 1, 2, 3, 4. c) Rescaled droplet profiles. The parameters 
of the simulation are At = 0.01, Ax = 0.05, b = 0.01. 



the fluctuations-dominated asymptotic rate of spreading. 

The appearance of a new power law for spreading sug- 
gests that the dynamics of the spreading drop is self- 
similar. The self-similarity is demonstrated through the 
excellent data collapse shown in Figs. 3b, 3c. Indeed, 
the power law £ s toch ~ t 1 ^ 4 follows directly by assum- 
ing a self-similar solution (in the statistical sense) of the 
form © to Eq. © [H- Conservation of drop volume 
in one-dimensional geometry implies f3 = 1 in the self- 
similar form ©■ Substituting this in Eq. © gives rise 
to two possible scalings. The first, obtained by assuming 
that the surface tension term is dominant on the RHS 
of Eq. ©, is just Tanner's law £det ~ t 1 ' 7 . By con- 
trast, if the stochastic force is dominant we consider a 
self-similar dynamics of the average (h) = f(\£\/t z )/\x\. 
On the RHS of Eq. © we substitute the average stochas- 
tic force (£) = 1/ \Jt\x\ over a time interval i and space 
interval Thus we obtain £ s t oc h ~ t 1 / 4 , in agreement 
with the asymptotic behavior of our simulation. Assum- 
ing the scaling relations h ~ and t ~ |i| 4 we eval- 
uate the surface tension term in Eq. © as \x\~ s and 
the average of the stochastic force as V2a\x\~ 5 . Re- 
quiring dominance of the stochastic term and returning 
to dimensional variables, we obtain that the stochastic 
scaling behavior is expected in the regime 



\x\ > x* 



h < hl/x* 



(11) 



where x* = h 7 /6 W 1/6 /e T 1/3 and £ T = y/k B T/j. Thus, 
a necessary condition for observation of stochastic scal- 
ing behavior is h ^> W ^> It- Typical values of £t are 
few angstroms for molecular fluids (far from the critical 
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point), or colloid size for colloidal suspensions PJ. Ap- 
plying similar analysis for the radially symmetric spread- 
ing, Eq. (JBJ, shows that in this case self-similarity dom- 
inated by the stochastic force gives rise to an enhanced 
rate of spreading £ s toch ~ t 1 ^ , compared to Tanner's law 
idet ~ t 1 / 10 . The average stochastic force is dominant 
over surface tension if 

r»^ /3 /V /3 , h«h 2 /3 £ T ^ , (12) 

and thus requires ho 3> It- 

For molecular fluids, the enhanced rates associated 
with stochastic scaling behavior can be observed if the 
stochastic force in Eq. @ is dominant not only with re- 
spect to the surface tension term as is expressed in Eqs. 
(I11I12|I , but also with respect to the VdW force that can 
be quite strong for thin films 13]. In the nondimensional 

units of Eq. @ this force takes the form ^7^(^- _1 f|)- 
Following a similar approach to the one that led to Eq. 
(|llfl we use for one-dimensional geometry the scaling re- 
lations h ~ and t ~ |ai| 4 to compare between the 
average stochastic and VdW forces and obtain the addi- 
tional condition for the stochastic scaling regime 

x « /i 3 / V /3 /W^ 2/3 ^ 1/6 , (13) 

where tvdw = yAjl- An overlap between the intervals 
in Eq. JHJ and (JT3J) is achieved if ho < W {£ v dw / ir) 2 ■ 
Typical values of A are 100/cbT p| and thus iydw > @t, 
and this overlap can be obtained for h ^> W ^> It- 
Similar analysis for two-dimensional geometry yields the 
result r <C {h^ir 2 / ivdw 4 ) 1 ^ 9 ■ Consistency of this con- 
dition with Eq. ljT3|) is possible only if ho <C ir 5 ' /ivdw , 
which seems unfeasible for typical fluids. We conclude 
that volume preserving fluctuations-dominated spreading 
can be observed for molecular fluids in one-dimensional 
flows if the initial height of the drop is large enough. 

By contrast, for complex fluids VdW forces with 
the solid plate are not expected to significantly affect 
the spreading process, and thus we require x (r) <C 
hg/a ,h 3> a for one- and two-dimensional flows, re- 
spectively. An overlap with Eqs. I|11I12|) is achieved 

if ho > (a 6 W/i T 2 ) 1/5 , a?/ 2 /t T 1/2 for one- and two- 
dimensional flow geometries, respectively. Both condi- 
tions are easily realized if the initial drop is large enough. 



Another spreading dynamics in which a stochastic scal- 
ing behavior might be observed is a "leaking" process, 
in which the height of the film at x = is fixed to a 
constant value ho by a continuous supply of fluid. A 
self-similar dynamics in this case has the form Q with 
(3 = 0. Following similar analysis we obtain for one- 
dimensional geometry the scaling behaviors: idet ~ t 1 ^ 4 , 
f-stoch ~ t 1 / 3 , while for two-dimensional geometry both 
idet and Istoch ~ t 1 / 4 . For the one-dimensional geom- 
etry we obtain an increased asymptotic rate of spread- 
ing due to thermal fluctuations. To check the realiz- 
ability of the stochastic scaling regime we compare the 
average stochastic force with surface tension and VdW 
terms, using the scaling relations: h ~ const and f ~ x 3 . 
For molecular fluids where VdW forces are important 
the stochastic scaling behavior is expected in the regime 
h^ /2 W^ 2 /£ T < x < hl /2 £ T /W 1 ' 2 i V dw 2 , which is again 
possible provided ho <C W(£vdw I^t) 2 ■ Stochastic scal- 
ing behavior of a complex fluid drop under this condition 
is achieved if h > a and x > h^^W 1 ^ 2 fix- 
To conclude, we derived a Langevin lubrication equa- 
tion @, and showed that it gives rise asymptotically 
to significant deviations from Tanner's law of spreading. 
By comparing the average stochastic force with classical 
forces such as surface tension and Van der Waals, we 
showed that fluctuation-assisted spreading is expected 
asymptotically in various flow geometries of molecular 
and complex fluids. Complex fluids are attractive candi- 
dates for studying this phenomenon, since confocal mi- 
croscopy techniques enable direct observation of macro- 
scopic thermal effects in systems such as colloidal sus- 
pensions. Direct imaging of thermal capillary waves of 
an interface between two colloidal liquids was recently 
achieved by such methods [l4|, and we believe that sim- 
ilar tools may enable observation of the enhanced rate 
of spreading predicted in this letter. We hope that this 
result, together with other examples, will motivate fur- 
ther studies of the role of thermal fluctuations in small 
dimensional fluid systems. 
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